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Josephson tunneling between YBa2Cua07_5 and Pb with the current flowing along the c-axis of 
the YBa2Cu307-a is presumed to come from an s-wave component of the superconductivity in 
YBa2Cu307„4. Experiments on multi-twin samples are not entirely consistent with this hypothesis. 
The sign changes of the s-wave order parameter across the Nt twin boundaries should give cancel- 
lations, resulting in a small (\?Nt) tunneling current. The actual current is larger than this. We 
present a theory of this unexpectedly large current based upon a surface effect: disorder-induced 
supression of the d-wave component at the (001) surface leads to s-wave coherence across the twin 
boundaries and a non-random tunneling current. We solve the case of an ordered array of d+s and 
d-s twins, and estimate that the twin size at which s-wave surface coherence occurs is consistent 
with typical sizes observed in experiments. In this picture, there is a phase difference of rr/2 between 
different surfaces of the material. We propose a corner junction experiment to test this picture. 

PACS Nos. 74.50. +r, 74.72. Bk 



I. INTRODUCTION 



Understanding the nature of the order parameter is 
one of the main challenges in the theory of high-T c su- 
perconductivity. One of the most fundamental issues 
is the competition between s-wave and d-wave. d-wave 
is definitely the rule for high-T c systems, yet there is 
strong evidence that the electron-doped NdCeCuO ma- 
terial is s-wave. At a microscopic level, this suggests 
that the mean-field pairing interaction has two eigenval- 
ues which vie for dominance. Understanding this compe- 
tition would provide insight into the pairing mechanism. 

In this regard YE^CusOy-a (YBCO) is of particular 
interest. In other systems, the square symmetry forces 
the ordering to be pure d-wave or pure s-wave. The 
presence of one most likely prevents the emergence of the 
other because of repulsive terms in the free energy, and 
the competition has a clear winner. In contrast, YBCO 
has orthorhombic symmetry. This makes it inevitable 
that d-wave and s-wave are always mixed 0. Joseph- 
son tunneling experiments with current flowing mainly 
in the a-b plane [^) have made it clear that the dominant 
component is d-wave, but a substantial body of work has 
also demonstrated the existence of Josephson tunneling 
along the c-axis from YBCO to a Pb electrode, indicative 
of an s-wave component || [jfj || Q. It is to be hoped 
that these latter experiments, if carefully analyzed, can 
tell us about the strength of the s-wave admixture in the 
order parameter. 

We shall first review the experimental and theoreti- 
cal situation concerning c-axis tunneling, concluding that 
there are major theoretical puzzles still to be resolved. 
Then we shall present a new model of the phenomena 
which we argue is in agreement with the data as they 
stand, and how to test the model more thoroughly. 



c-axis tunneling from YBCO to Pb was first observed 
in twinned crystals etched wth Br ||, with an I c R n 
product of as much as 10 % of the known gap of about 
30 meV. This strongly suggested the presence of an s- 
wave component of the superconductivity of YBCO, as 
a pure d-wave current would average to zero over the 
Fermi surface. However, another possibility was that the 
current was due to second order tunneling of the d-wave 
component 0. This hypothesis predicts the presence 
of Shapiro steps in the conductivity in units of /i//4e, 
where / is the frequency of the incident radiation. This 
was ruled out in subsequent microwave experiments Q 
Finally, the question of tunneling through step walls at 
the surface arises, particularly if it is deeply etched. This 
would be a process in which the current actually flows in 
the a — b direction. However, the fabrication of a — 
b junctions Q and the observation of tunneling in situ 
without etching [^| appears to have laid this possibility 
to rest. The presence of a nonzero s-wave component in 
YBCO must now be accepted. 

Is it reasonable to accept the 10% estimate of s to 
d which comes from the I c R n product at face value? 
Clearly not, for the following reason. A twinned sam- 
ple should have a relative population of twins of the two 
possible orientations given by statistical considerations. 
The d-wave component must remain coherent through 
the sample, as is shown by the corner junction experi- 
ments B. Because the change in orientation reverses 
the relative sign of s and d we should have roughly equal 
areas of d + s and d — s superconductivity in the sample, 
in which case the net current should be zero. More 
precisely, the net current should be proportional to to 
I c \/Nt, where I c is the critical current of a single twin 
and Nt is the total number of twins. If we accept this 
argument, then the actual proportion of s to d would be 
higher than 10%. This would move the nodes in the 
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gap well away from the diagonal in the Brillouin zone. 
This would be inconsistent with tricrystal experiments 
Pj Furthermore, comparison of Josephson currents in 
single crystals to twinned crystals show similar R n values 
and I c values which range from 0.5 to 1.6 mA for single 
crystals and from 0.1 to 0.9 mA for twinned samples [Q. 
These numbers are subject to the objection that one can- 
not be sure that the tunneling matrix elements are not 
extremely sensitive to the sample preparation method. 
Nevertheless, in view of the fact that R n does not vary 
wildly from sample to sample, they suggest that a purely 
statistical analysis of twin populations with a resulting 
small imbalance of d + s and d — s is not a viable expla- 
nation of the data. 

The dilemma was deepened by experiments on crys- 
tals with much larger twins, large enough so that junc- 
tions could be formed which straddled either one or even 
zero twin boundaries These showed that the direc- 
tion of current definitely did change sign across the twin 
boundary, a fact which can be established unambiguously 
by investigating the current as a function of field in the 
plane of the junction. This observation was consistent 
in all eight samples studied. Also, no such sign changes 
were observed in the absence of a boundary. These ex- 
periments therefore clearly confirm the mechanism of an 
s-wave component controlled by the orthorhombic distor- 
tion, without offering any explanation of the large current 
in heavily twinned samples. One further observation in 
these experiments may offer a clue, however. The cur- 
rent at zero appplied field in single-boundary samples 
was consistently higher than calculated by looking at the 
relative sizes of the two twins. We will return to this 
point below. 

Summarizing the experiments, we may say that an s- 
wave component which changes sign across boundaries is 
clearly present. If it always changes sign, then we can- 
not explain the data on twinned samples using purely 
statistical arguments. One possibility is that the twin 
populations are not equally likely. For example, if the 
twinning takes place under uniaxial stress, then one ori- 
entation would be favored. Experiments which correlate 
microstructure with Josephson current are needed to rule 
this out Q|. However, given the size of the Josephson 
effect in twinned samples, it seems to us that this expla- 
nation is somewhat implausible. 

The most detailed theory of c-axis tunneling proposed 
to date is that offered by Sigrist et al. . Their picture 
involves no net tunneling from the twins themselves. A 
time-reversal-breaking state at the twin boundary is pre- 
dicted which results in a net Josephson current coming 
from the twin boundaries. This would give a Josephson 
current which is proportional to the number of bound- 
aries for a fixed surface area. This is not observed, 
though again one must keep in mind that different sam- 
ples must be compared to make any such statement, and 
variations in important microscopic parameters canot be 
controlled in such comparisons. In addition, however, 
the theory predicts a current which has maximum asym- 



metry (as a function of in-plane angle) when the applied 
magnetic field is parallel to the boundary. This is an 
experiment in which the unknown matrix elements are 
held fixed. This prediction is in conflict with the ex- 
perimental observations, which are symmetric at parallel 
orientation |J. 

A quite different proposal was made by Xu et al. pl| . 
These authors postulate a bulk d + is state. In this 
theory, however, the s-component does not change sign 
across the boundary, which does not agree with the mea- 
surements on single-boundary samples in a parallel field. 

We present an alternative explanation in which the 
nonzero tunneling current is the result of a surface ef- 
fect. YBCO is notorious in photoemission experiments 
for not showing a gap. This proposal is inspired by the 
fact that photoemission experiments (with resolutions of 
order WmeV or less) have also never succeeded in seeing 
a gap at the (001) surface in this material (in contrast 
to Bi 2 Sr2CaCu 2 08+x)- This shows that the magnitude 
of the gap at this particular surface is much reduced. 
Furthermore, if this reduction is due to disorder, such 
as surface scattering, one would expect that the d-wave 
component is relatively much more suppressed than any 
s-wave admixture. A similar suppression could result 
from an oxygen vacancy concentration gradient. This 
suppression of the d-wave component of the order pa- 
rameter as we approach the (001) surface of the YBCO 
in the context of c-axis tunneling is of two central hy- 
potheses of our model and was first suggested by Bahcall 
fll2fl . The second crucial ingredient is new: the d-wave 
surface suppression results in a coherent s-wave surface 
layer and hence an enhanced Josephson tunneling current 
in highly twinned samples without a very large admixture 
of s-wave in the bulk. 



II. DOUBLE TWIN MODEL 

Twinned samples are disordered on the /im scale, the 
twin boundaries running predominantly along the diag- 
onal of the a-b plane. It is reasonable then to approxi- 
mate the disordered sample by an array of straight twin 
boundaries running across the entire sample, which is 
considered to be semi-infinite. We concern ourselves in 
this paper with the ordered case in which all twins are of 
the same width and alternate between d+s and d-s. In 
real samples the twins have varying widths, but we have 
verified numerically that the basic results are unaffected 
by neglecting the disorder in the widths. The solution 
of the ordered model should be periodic with a period 
of two twins. Therefore, we solve the case of two twins 
with periodic boundary conditions. The twin boundary 
occupies the half plane defined by x — and z < 0. The 
plane z = is the (001) surface of the YBCO sample. 
The model is illustrated schematically in Fig. [l]. 
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FIG. 1. The ordered case where all twins are of equal width 
is equivalent to a d+s and a d-s twin with wrap around bound- 
ary conditions. Near the (001) surface, there is a layer in 
which the d component of the order parameter is highly sur- 
pressed. 



We shall write the bulk free energy density in terms 
of the two order parameters ^ d and ty s . Spatial varia- 
tion is allowed only along the x-direction, (normal to the 
twin boundary), and along the z-direction (normal to the 
surface). The free energy density is given by: 



f = a d (z)\y c 
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V d \ 4 + K dx \d x * d \ 2 +K dz \d z y d \ 



+a s \V s \ 2 + ^|* s | 4 + K sx \d x * s \ 2 + K SZ \8 Z * S \ 
+a sd (x){^ d + c.c) + l3 sd \y s \ 2 \V d \ 2 



(1) 



Some of the important physical ideas behind our model 
are displayed by this equation. a d is a function of 
position in order to enforce the condition that the d- 
wave component is suppressed near the surface. Thus 
o-d — * &do < 0, a negative constant, as z — > — oo, deep in 
the bulk. a d increases toward as z — > at the surface. 
a sd is the s-d coupling parameter which is a negative con- 
stant deep in the d + s twin and positive constant deep 
in the d — s twin. a sd {x) = —a sd (—x). 

The (3 sd term is the s-d repulsion mentioned in the 
introduction. We will neglect it in the calculations and 
have included it here only in order to stress that a large 
positive (3 sd suppresses all s-d mixing in the absence of 
the bilinear a sd term. If this term is present, as it is here 
because of the orthorhombic distortion, then the size of 
the s admixture is controlled by a sd /a d . 

We must also include the free energy of the twin bound- 
aries. Any x axis variation of \I/ S and ^ d will take place 
within a distance of the order of a coherence length about 
the twin boundary. Since £ a b ~ 20A is very small com- 
pared to the average twin width (0.1 to 10/ito) we con- 
clude that the detailed structure of the twin boundary is 
not very important. We will assume a very thin bound- 
ary and thus take a sd to be piecewise constant. This is in 
direct contrast to the Sigrist et al. model in which the cur- 
rent comes from the twin boundaries. In our model the 



current comes from the twins. We therefore approximate 
the free energy of the twin boundary by a Josephson-type 
coupling: 
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where J s and J d > 0. We can also now drop the x axis 
gradient terms in the bulk free energy. Any x axis vari- 
ation in the order parameters takes place near the twin 
boundary and has been included in the boundary energy. 

The problem has been reduced to 2 one-dimensional 
twins which are Josephson coupled. However, only one 
twin is actually required. As the surface of the YBCO 
is approached, the magnitude of ty d and \P S should vary 
in exactly the same way in both the d + s and d — s 
twins. Only the phases <j> s and <j) d arc different. But while 
the phases differ between twins, they are not entirely 
independent. We set <p d = 4> s = in the bulk of the d + s 
twin, and (f) d = 0, cf) s — ir in the bulk of the d—s twin. As 
the (001) surface is approached, variation in </>+ and <j>j 
should be symmetric about n/2, while <j)^ and (f)^ will be 
symmetric about 0. This allows the boundary energy to 
be rewritten entirely in terms of the phases in the d+s 
twin. 



rt+ s (z) = -J d \* d (z)\ 2 cos(2Mz)) 
+ J s \^ s (z)\ 2 cos(2</> s (z)) 



(3) 



<f> s and <f> d in the d—s twin can be deduced immediately, 
and the problem is now entirely one-dimensional. 



III. SOLUTION 

We will solve for the order parameters in the d+s twin. 
The solution for the d—s twin follows immediately. Our 
one-dimensional free energy density is: 

f(z) = w{a d (z)\* d \ 2 + ^ d \ 4 + a s \* s \ 2 + ^|* s | 4 

+ K dz \d z * d \ 2 +K sz \d z * s \ 2 
+ a sd {x)(^*^ d + c.c)} 

+ -J d |* d | 2 cos(20 d ) + J s |* s | 2 cos(20 s ) (4) 

where w is the width of a single twin. Performing the 
usual minimizations, we get: 

Sf 
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w{a d (z)\V d \ + /3 d \V d \ 3 + a sd \^ s \ cos(0 d - <f> s ) 



and 
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25<f> d 



+ K dz (--d 2 z \y d \ + \y d \ 2 (d z M 2 )} 

- J d |* d |cos(20 d ) = 



w{-a sd \^ d \\^ s \sm((f) d - cj) s ) 



(5) 



-K dz \^ d \ 2 d 2 <t> d ] + J d \^ d \ 2 sin(20 d ) = (6) 
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The analogous s-wave equations are: 
Sf 



25\t> s 



and 



Sf 
2Stp s 



= w{a s \^ s \ +(3 s \^ s [ i + a sd \^d\cos( 

+ K sz (-^d 2 \^ s \ + \^ s \ 2 (d z cj )s ) 2 )} 
+ J s |# s |cos(2<£ s ) = 



= iu{-a s d|*<i||*s| sin(<?!> s - cj> d ) 



(7) 



-K sz \<£ s \ 2 d 2 ^ s } - J s |^ s | 2 sin(20 s ) = (8) 



In general, these equations must be solved numer- 
ically, but it is instructive to first consider the limit 
Kdz = Kds = 0, which may be obtained analytically. 
If we consider Eqns. 6 and 8 we see that: 



J d \^ d \ 2 sin(20 d ) = J s |* s | 2 sin(20 5 .) 
= iua s d|* s ||*<i| sin(0 s - <j> d ) 



(9) 



(j> s and <f>d are between and 7r/2. There are only the 
two obvious solutions: <fid = <j> s = or tt/2. 

The particular solution which minimizes the free en- 
ergy is dependent upon the relative strengths of the s 
and d intertwin Josephson couplings, i.e. , the ratio 
R = Jd|*d| 2 M|* s | 2 . H R > 1 then (f> s = <j) d = 0. If 
R < 1 4> s — 4>d — tt/2. This is in the d + s twin. In the 
d—s twin if R > 1 we have 4>d = 4> s = tt. If R > then 
<f>d = —tt/2 and <j) s = tt/2. The magnitudes are obtained 
from the coupled set of equations: 
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\a sd \ 
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\a s d\ 



-{a d (z)\y d \+/3 d \V d \ 3 ) 
f (a s |^ s |+/? s |^ s | 3 ) 
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where we have assumed that the intertwin coupling has 
little effect on the magnitude of the order parameters, 
that is J s << w\a s \ where w is the twin width. The 
exact form of |^ s | and \^ d \ will depend on the a d (z) 
chosen. 

The main effect of finite gradient terms K s d \d z ^ s ,d\ 2 
in the free energy is to smooth out the variation in the 
order parameters as the surface is approached. We ex- 
pect the order parameter magnitudes to be only slightly 
affected by the introduction of the gradient terms. Vari- 
ation in d\ and |\I/ S | should depend predominantly on 
a d (z), since K d « \a d \ etc. The effect on the phases 
is more dramatic. For relatively narrow twins, <f> s and 
(j)d now undergo a smooth transition from <p s — <p d = 
in the bulk of the d + s twin to <f> s = (f> d = tt/2 at the 
surface. In the d—s twin, <j) s changes from tt to tt/2 at 
the surface and <j) d from to —tt/2. The order parameter 
magnitudes and phases for a model a d {z) are shown in 

Fig. |. 
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FIG. 2. Order parameter magnitudes and phases near 
the (001) surface (z=0). a d (z) = -0.95(1 - e z ) + 0.05, 
a 3 = -0.5, (i d = (3s = 0.2, K d = K a = 10- 6 , \a sd \ = 0.01, 
Jd — Js — 0.005. In the bulk <j>d is coherent across the twin 
boundary. At the surface, <f> s is coherent and is out of phase 
with the bulk 4>d by 7r/2. 



The degree of smoothing depends upon the strength of 
the gradient term versus that of the coupling across the 
twin boundary. The dominant factor in this competition 
between the gradient and the intertwin coupling energy 
is the twin width. For very wide twins, the change in 
surface phase is diminished and may be eliminated alto- 
gether. 

The maximum c-axis Josephson current is given by: 
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(11) 



where A is the junction area and (j)pb has been cho- 
sen to yield the maximum Josephson current. For very 
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large twins it is not energetically favorable for the phase 
change to occur. The s-wave phase at the surface al- 
ternates between and it across twin boundaries and 
no net Josephson current flows. As the twins become 
narrower, a threshold is reached where the s-wave phases 
start to shift towards ir/2 at the surface. The s-wave sur- 
face phase alternates between </>f +s and </>f~ s = tt — <j>s +s - 
Some Josephson coupling is now possible. For very nar- 
row twins 4>s is coherent across the entire surface of the 
crystal, and the maximum Josephson current flows. The 
current saturates, and further reduction of twin size has 
no effect on the current. This is illustrated in Fig. ||. 
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FIG. 3. Net Josephson current vs twin size. 
Landau parameters are as in Fig 2. J = J s 
K = K a = Kd- T is the twin size. 
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The saturation is one important phenomenological dif- 
ference between our model and that of Sigrist et ai. 

We want an estimate of the average twin width at 
which s wave surface coherence begins in terms of ex- 
perimentally measurable quantities. Roughly speaking, 
this threshold will occur when the strength of the s-wave 
coupling between twins is equal to the gradient energy in- 
volved in rotating the s-wave phase by tt/2 at the surface. 
We will assume an very simple model with a surface layer 
of depth s in which |* d | = |\f^| for z < —s and |* d | = 
for — s < z < 0. \^ s \ is assumed constant for all z. 

The first task is to get some idea of the strength of 
s-wave coupling across the twin boundary. If we consider 
the situation far from the surface we may take 1^(2 = 
— oo)| to be large and fixed. An effective free energy 
may then be written down for ty s and an Euler-Lagrange 
equation for the variation in *5f s with respect to x derived. 



a s * s + a sd ^ d - K s 



d 2 ^ s 
dx 2 



= 



(12) 



If we assume a step function boundary where a s d(x) 
—sgn(x)a sd then we have the solution 



9 B (x) = sgn(x)^° s (l - e~^) 



(13) 



where ^ is the bulk value. The result for the free energy 
per unit area of the twin boundary is: 



Fb 
A 



I 2 



(14) 



The c-axis gradient energy is also required and is 
roughly 



5L=K s \* s \->C-llf W 



(15) 



where w is the twin width. Noting that K s /a s = and 
setting F g = F b we obtain: 



Sab 2 
Sc 



(16) 



For a surface layer with a depth of lOOA (about 8 unit 
cells) then we obtain a twin width of approximately ljum. 
We emphasize that this is merely an order of magnitude 
estimate. In addition, it is not clear exactly how deep 
such a surface layer should be. However, the resulting 
twin width is not unreasonable. A highly twinned sample 
of 0.5mm may have as more than 10 3 twins resulting in an 
average twin width of a few tenths of microns. Thus while 
we expect no net Josephson current in a lightly twinned 
sample, our model predicts the net Josphson current ob- 
served in more heavily twinned crystals. 



IV. PROPOSED EXPERIMENT 

Our model predicts a nonzero Josephson current re- 
sulting from a surface effect. For samples with relatively 
large twins, we expect a d+s d-s alternation between 
twins at the surface and no net Josephson current. This 
explains why experiments on two twin crystals show a 
sign change in the Josephson coupling to Pb across the 
twin boundary 0]. In a sample with many smaller twins, 
however, the coupling between twins wins out and a co- 
herent s-wave surface layer results. We expect this to 
take place in samples where the average twin width is 
less than a few micrometers. The s-wave surface layer is 
7r/2 out of phase with the bulk d-wave phase. 

We emphasize this fact because the 7r/2 phase shift is 
experimentally verifiable. A YBCO-Pb corner junction 
type experiment with one junction on the (100) surface 
and the other on the (001) surface of a highly twinned 
YBCO sample should be able to detect this 7r/2 phase 
shift, as was previously suggested by Sigrist et al. |jTo| . 
We give a schematic diagram of the proposed experimen- 
tal configuration in Fig. ^. The current maximum as a 
function of field will be shifted by a quarter of a flux 
quantum. The Josephson coupling to the Pb at the (100) 
junction is predominantly due to the YBCO d-wave com- 
ponent since d-wave suppression is not expected at this 
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surface. Since the c-axis Josephson coupling results from 
the smaller s-wave component, it is much weaker than 
the a-axis coupling. The (001) junction should therefore 
have a much larger area than the (100) junction in order 
to minimize any DC offset of the interference pattern. 
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FIG. 4. Schematic of proposed YBCO-Pb SQUID experi- 
ment. One junction is on the (100) surface of the YBCO, the 
other on the (001) surface. The tv/2 shift between the bulk 
d-wave phase and the (001) surface s-wave phase shifts the 
current maximum by a quarter of a flux quantum. 



V. CONCLUSION 

The present theory can reconcile the puzzles mentioned 
at the outset. The suprprisingly large value of I c R n is as- 
cribed to the partial, coherence of the s-wave component 
at the surface. The fact this coherence is only partial 
gives a reasonable account of the overall differences be- 
tween single crystal and twinned samples. The fact that 
single-boundary junctions always show a change in sign 
of the s component is also consistent: in this case the 
twins are larger. These experiments also show excess 
current at zero applied magnetic field. This would be 
consistent with some partial coherence of the s compo- 
nent across the boundary, as the larger (stronger) of the 
two twins appears to control the weaker one. Hence we 



believe that the theory can account for all observations. 
The experiment of the previous section would be a criti- 
cal test of the theory. Experiments in which the relative 
twin populations are precisely controlled would serve to 
rule out the alternative explanation in which the current 
is due to accidental anisotropy introduced in the growth 
process. 

One important qualitative conclusion about the under- 
lying physics of the bulk can be drawn from this picture: 
s-wave competes with d-wave in YBCO. If our model 
is correct, then the naive estimate of 10% admixture of 
s-wave as a proportion of d-wave remains roughly cor- 
rect. Expressed in the language of Eq. [jj we have that 
|*s|/|*d| °" a sc i/{a s — ad) °" 0.1 at low temperatures. 
If s-wave were very strongly suppressed by a large posi- 
tive a s , it would not be so easily induced by the lattice 
distortion. 

The present theory would predict that only those ma- 
terials with orthorhombic distortion should show c-axis 
tunneling. Recently, c-axis Josephson tunnneling be- 
tween Ba2Sr2CaCu208+a; and Pb [|13| has been observed 
in spite of the absence of an orthorhombic distortion in 
this material. However, due to the fact that I c R n ~ 
lfieV, orders of magnitude less than the gap value, we 
believe that this interesting effect is physically different 
from that seen in the YBCO experiment. 

We would like to thank J. Betouras for many helpful 
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under the Materials Research Science and Engineering 
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